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ABSTRACT 

We propose a scenario for the formation of DA white dwarfs with very thin helium 
buffers. For these stars we explore the possible occurrence of diffusion-induced CNO- 
flashes, during their early cooling stage. In order to obtain very thin helium buffers, 
we simulate the formation of low mass remnants through an AGB final/late thermal 
pulse (AFTP/LTP scenario). Then we calculate the consequent white dwarf cooling 
evolution by means of a consistent treatment of element diffusion and nuclear burning. 

Based on physically sounding white dwarf models, we find that the range of helium 
buffer masses for these diffusion-induced novas to occur is significantly smaller than 
that predicted by the only previous study of this scenario. As a matter of fact, we find 
that these flashes do occur only in some low-mass (M < 0.6A/ Q ) and low metallicity 

(Zzams <; 0.001) remnants about 10 6 — I0 7 yr after departing from the AGB. For 
these objects, we expect the luminosity to increase by about 4 orders of magnitude in 
less than a decade. We also show that diffusion-induced novas should display a very 
typical eruption lightcurve, with an increase of about a few magnitudes per year before 
reaching a maximum of My ~ -5 to -6. Our simulations show that surface abundances 
after the outburst are characterized by logA^H/-/VHc ~ —0.15. ..0.6 and N>C;>0 by 
mass fractions. Contrary to previous speculations we show that these events are not 
recurrent and do not change substantially the final H-content of the cool (DA) white 
dwarf. 

Finally, with the aid of model predictions we discuss the possibility that Nova 
Vul 1670 (CK Vul) and the recently proposed [WN/WC]-central stars of planetary 
nebulae could be observational counterparts of this diffusion-induced nova scenario. 
We conclude that, despite discrepancies with observations, the scenario offers one of 
the best available explanations for CK Vul and, with minor modifications, explains 
the observed properties of [WN/WC] -central stars of planetary nebulae. 

Key words: stars: AGB and post-AGB -stars: individual: CK Vul, PB8 -stars: Wolf 
Rayet - white dwarfs. 



1 INTRODUCTION 

It is well known that element diffusion plays a major 
role in shaping the chemical stratification of white dwarfs 
l|Koester fc Chanmu gam 1990). Gravitational settling is re- 
sponsible for the chemical purity of the outer layers of white 
dwarfs. Due to the very high surface gravities of these stars 
the lightest element remaining in the white dwarf (either H 
or He) floats up and forms a pure envelope, the thickness 
of which gradually increases as the white dwarf evolves. On 
the other hand, the large chemical gradients formed dur- 
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ing the previous evolution, lead to chemical diffusion which, 
in turn, wil l smooth out the ch emical profile of cool white 
dwarfs (e.g. lRenedo et alj|2010l) . Recently, the gravitational 
settling of 22 Ne in the liquid interior of white dwarfs, has 
been found to be a key ingredient to explain t he cooling age 
of the old, metal-rich, open clus ter NGC 6791 (jAlthaus et al.1 
l2010l . iGarci'a-Berro et al.ll201Gn . 

Since the first simulations of white dwarf evolution that 
incl uded a simultaneous treatment of diffusion and cooling 
(i.e. Ilben fc MacDonaldlll985l . [ 1986) it was noticed that dif- 
fusion could trigger thermonuclear CNO-flashes. In fact, the 
inward diffusion of H and the outward diffusion of C within 
the pure He zone that is located below the H-rich envelope 
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and above the C-rich intershell (usually named "He- buffer" , 
see Fig. [T]) can lead t o a runaway CNO-burnin g. This pos- 
sibility was shown bv llben fc MacDonaldl l|l986h who found 
that, if the He-buffer was thin enough, C and H can come 
into contact and develop an unstable H-burning. This sud- 
den energy release produces a very rapid expansion, of the 
order of years, of the outer layers of the white dwarf push- 
ing the star back to a giant configuration and increasing its 
visual magnitude from My ~ 9 to My ~ —6 in a few years. 
We term this eruptive event as "diffusion-induced nova'f]] 
(DIN) although it leads to a much s lower brightening than 
classical novas. Ilben fc MacDonaldl (|l986h showed that af- 
ter such events the stars will become, in a few years, yellow 
giants with mildly He-enriched surface compositions. In a 
more speculative mood they also suggested that DINs may 
be recurrent, finally leading to H- deficient compositions. 
Later, prompted by this speculation. iD'Antona fc Mazzitellil 
( 1990) suggested that the H-rich envelope could be strongly 
reduced during these events, leading to DA white dwarfs 
with thin H-envelopes, as inf erred in some DA white dwarfs 
i|Castanheira fc Keplerll2009l ). 

Since the first explo ratory simulation performed by 
Ilben fc MacDonaldl (1986), these DINs have been mentioned 
in sev e ral review p a pers and book s llD'Antona fc Mazzitelli 
1 1991 llbenl 1 19951 . 
and more recently 
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ISalaris fc Cassisil 12005 
|2010| ). but no new nu- 
merical experiments have been performed. The latter 
becomes particularly relevant i n ligh t of the criticism 
made by iMazzitelli fc D'Antona ll 19861') that the pre- white 
dwarf models adopted by Iben fc MacDonaldl |l985l ) and 
Ilben fc MacDon ald (1986) were not realistic enough. In this 
context, a reexamination of the possibility and characteris- 
tics of DINs is in order. 

As mentioned, DINs are expected to increase their vi- 
sual magnitudes from My ~ 9 to My < —5 in a few years. 
Ilben fc MacDonaldl l| 19861 ) mentioned that such event would 
display a lightcurve similar to that of a very slow nova. The 
most famous object with such description is the enigmatic 
Nova Vul 1670 (CK Vul), the oldest catalogued nova variable 
|Trimblell2000l ) . First discovered by Pere Dom Anthelme the 
20th of June of 1670, and later independently by Johannes 
Hevelius on the 25th of July, seem to have rai sed to 3rd 
magnitude in less than a year (|Shara et al.lll985l l. although 
no pre- maximum observations exist. Its eruption lasted for 
two years (1670-1672). During that period it faded and 
brightened two times befo re definitely disappea ring from 
view. Three centuries later. IShara fc Moffat! (|1982T ) partially 
recovered the remnant. They discovered several nebulosi- 
ties that were concl usively linked to the 1670 eruption by 
lHaiduk et ail ([20071 ) through proper motion studies. CK 
Vul remains a very enigmatic object and several scenar- 
ios have been proposed to explain its observed behaviour. 
In fact, it has alternatively proposed to be an hibernating 
nova, a very late thermal pulse (VLTP), a " gentle" super- 
nova, a light nova and a stellar merger event (|Haiduk et al.l 



Ilben fc MacDonal d (1986) named these events self-induced no- 
vas, but we will use the term "diffusion-induced novas" to distin- 
guish them from late helium flashes (born again AGB stars, llbenl 
1 19841) which predict similar lightcurves and can also be considered 
"self-induced" . 



l2007[|Katcll2003l ). None of the pro posed scenario seem s to 
correctly describe the observations (jHaiduk et al.ll2~007l). al- 
though the VLTP seem s t o be the preferred o ne ( Harrison! 
1 19961. IE vans et al.ll2002l and lHaiduk et al.ll2007f ). In this con- 
text, and due to the similarities between post- VLTP evo- 
lution and post-DIN evolution, it seems worth analysing 
whether CK Vul could be understood within the DIN sce- 
nario. 

The main purpose of the present article is to study the 
possibility that DINs could take place in physically sound- 
ing white dwarf models with a realistic evolutionary history. 
We will also identify a detailed scenario for the creation of 
white dwarfs with thin enough He-buffers for DIN events 
to occur. Specifically, to perform this study we compute re- 
alistic white dwarf models by means of "cradle to grave" 
stellar evolution simulations. Then we compute white dwarf 
cooling sequences by considering a simultaneous treatment 
of element diffusion and evolution. 

The article is organized as follows. In the next section 
we briefly describe the input physics adopted in the present 
simulations. Then, in section 3, we identify a possible evo- 
lutionary scenario that can lead to white dwarfs with thin 
He-buffers. Once the evolutionary scenario has been set we 
present our numerical simulations and results (Section 4). 
In particular, we show that the proposed scenario can ef- 
fectively produce white dwarf models with thin He-buffers. 
Then, we discuss the main predictions of the simulations and 
explore under which conditions DINs might occur. Further- 
more, we make a very detailed description of the sequence 
of events before, during and after the DIN events. Once 
the predictions of DIN simulations have been clarified, in 
Section 5 we look for possible observational counterparts of 
these events. In particular, we speculate that DINs could be 
the origin of CK Vul and PB8 (t he prototype of th e recently 
proposed [WN/WC]-CSPNe, see lTodt et atll2010h and show 
that, although significant discrepancies exist, the DIN sce- 
nario offers a possible explanation for these objects. Finally, 
we close the article with some final statements and conclu- 



2 BRIEF DESCRIPTION OF THE STELLAR 
EVOLUTION CODE 

The simulations presented here have been performed wit h 
LPC0DE stellar evolution code (see iRenedo et al. I l2010l ). 
LPC0DE is a well tested and amply used evolutionary code 
that is particularly appropriate to compute realistic white 
dwarf structures and their evolution by means of "cradle to 
grave" stellar evolution simulations, and allows for a con- 
sistent and detailed computation of the effects of element 
diffusion during the white dwarf cooling phase. As this work 
is aimed at exploring the outburst properties of diffusion- 
induced CNO-flashes on hot white dwarfs, it is worth men- 
tioning that LPC0DE has been previously used to study other 
kind of H-flashes on (very) hot white-dwarfs. These in- 
clude the simulations of the post- A GB born again scenario 
l|Miller Bertolami fc Althausll2007h. and the post-RGB hot- 



flasher scenario l|Miller Bertolami et al. 2008) for the forma- 
tion of H-deficient white dwarfs and subdwarfs respectively. 
Numerical and phys ical aspects of LPCODE can be found in 
lAlthaus et alj l|2005h . In what follows we only comment on 
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Figure 1. Sketch of a Kippenhahn diagram of the proposed sce- 
nario for the formation of DA white dwarfs with thin He-buffers. 



the physical inputs regarding the treatment of diffusion pro- 
cesses that are specifically relevant for the present work. 

White dwarf models are obtained as the result of com- 
puting the evolution of low mass stars from the ZAMS 
through the helium core flash and through the thermal 
pulses on the AGB (TP-AGB) and, then, to the white dwarfs 
stage. LPCODE considers a simultaneous treatment of non- 
instantaneous mixing and burning of elements, by means of 
a diff usion picture of conve ction coupled to nuclear burning 
— see lAlthaus et all (|2005h for numerical procedures. The 
nuclear network considered in the present work accounts ex- 
plicitly for the following 16 elements: 1 H, 2 H, 3 He, 4 He, 7 Li, 
7 Be, 12 C, 13 C, 14 N, 15 N, 16 0, 17 0, 18 0, 19 F, 20 Ne and 22 Ne, 
together with 34 thermonuclear reaction rates correspond- 
ing to the pp-chains, t he CNO bi-cycle, He-bu r ning a nd C- 
ignition as described in iMiller Bertolami et all (|2006h . 

As mentioned in the introduction, element diffusion 
strongly modifies the chemical profiles throughout the outer 
layers of white dwarfs. In this work white dwarf evolution 
was computed in a consistent way with the changes of chemi- 
cal abundances as a consequence of diffusion. The tr eatment 
is similar to that of llben fc MacDonaldl (|l985l , ll98tj ) but we 
consider, in addition to gravitational settling and chemical 
diffusion, the process of thermal diffusion. We do not take 
into account radiative levitation, as it is only relevant for 
determining surface chemical abundances and, thus, is irrel- 
evant for the purpose of the present work. Our treatment 
of time-dependent element diffusion is ba sed on the multi- 
component gas picture of iBurgersl (1 19691 ). Specifically, we 
solved the diff usion equations wi thin the numerical scheme 
described in lAlthaus et ail l|2003h . 



3 THE EVOLUTIONARY SCENARIO 

The formation of a WD with a very thin He-buffer seems 
to be a key factor i n the o ccurence of a DIN e pisode (see 
llben fc MacDonaldl Il986t ). Ilben fc MacDonaldl (| 19861 ) ob- 
tained a pre-white dwarf model with a thin He-buffer by 
abstracting mass from an A GB star model shortly a fter it 
has experienced a He-flash. iMazzitelli fc D'Antonal (|l986l ) 



noted that the way this pre-white dwarf model was created 
is artificial and could affect the thermal properties of the 
remnant. Whether white dwarfs with thin He-buffer can be 
actually formed relies on identifying a scenario in which they 
could be formed under standard assumptions. As shown in 
Fig. \T\ after a He-shell flash on the TP-AGB, the mass of 
the He-buffer region becomes strongly reduced by intershell 
convection. Such thin He-buffer survives until the reignition 
of the H-burning shell. Note that this may not be true if 
the star experiences third dredge up events, during which 
the deepening of the H-rich envelope will erase the remain- 
ing He-buffer by diluting it into the envelope. From this 
moment on, the He-buffer mass starts to increase until the 
next He-shell flash. In this context, a realistic scenario that 
can lead to a DA white dwarf with a thin He-buffer can 
be outlined. Recent studies on the initi al-final-mass relation 
(|Kalirai et al.l l2008. Salar is et al.ll20Qgf ) support th e absense 
of sig nificant third dredge up in low-mass stars (|Karakasl 
120031 ). In particular, the first thermal pulses of low-mass 
stars (M <J 1.5Mq) are not very strong and, thus, no third 
dredge up takes place in numerical models. Hence, it is not 
unreasonable to accept that low-mass stars experiencing an 
AGB final thermal pulse (AFTP) or a Late Thermal Pulse 
(LTP) will end as DA white dwarfs with thin He-buffers. In 
those cases, as no third dredge up happens, the very thin 
He-buffer survives the last helium shell flash (either AFTP 
or LTP). Then, during the He-burning phase that follows 
the flash, AGB winds will erode an important fraction of 
the remaining (already depleted) H-rich envelope, prevent- 
ing a reignition of the H-burning shell and an increase in 
the He buffer mass. As a result, the He-buffer is still very 
thin when the star finally reaches the white dwarf phase. 
We propose that this scenario can lead to white dwarfs with 
thin He-buffers and therefore might lead to DINs. 

In order to bear out the scenario described in the 
preceeding paragraph we have computed several sequences 
from the Z AMS to the TP-AG B under standard assump- 
tions — see iRenedo et ail (|2010h for a description. Partic- 
ularly, we have included overshooting only in core con- 
vective zones within a diffusive convective picture as in 
IMiller Bertolami fc Althausl (2006). Thus, no overshooting 
mixing has been included during the TP-AGB evolution of 
any sequence, as it leads to strong third dredge up episodes 
in low-mass AGB stars. Although the efficiency of third 
dredge up on AGB stars in clearly not a solved problem, it 
is clear that third dredge up is less e fficient in low mass s tars 
than in the mor e massive stars (e.g. iKarakas et al.ll2002l ) . In 
this connection. ISalaris et al.l (|2009l ) have shown that strong 
third dredge up is disfavoured, by the semiempirical initial fi- 
nal mass relation, for stars with masses lower than 4Mq. It is 
worth noting that their study only deals with stars of masses 
M > 2Mq and can not be directly applied to our work. How- 
ever, as third dredge up efficiency is expected to decrease for 
stars of lower masses, it seems safe to assume that no third 
dredge up should be pr esent in stars of even low er masses. 
In agreement with this, iGuandalini et ail (|2006) have sug- 
gested that the properties of carbon stars can be explained if 
no overshooting is added in low ma ss AGB stella r mode ls. It 
must be noted, h o wever , that both lSalaris et all (120091 ) and 
IGuandalini et al] l|2006l ) deal with stars of higher metallici- 
ties than studied in the present work and third dredge up is 
expected to increase at lower metallicities. The absence of 
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Figure 2. HR-diagram evolution of an initially O.85M0 star with Z = 1 X 10 — 4 . It undergoes its final thermal pulse on the AGB with 
a thin envelope and consequently no important H-reignition takes place before reaching the white dwarf cooling stage. This produces a 
white dwarf with a very thin He-buffer which experiences a DIN due to the diffusive mixing of C and H in the He- buffer (sec Fig. [3J. 
Time has been choosen to be t = at the point of maximum energy release by H-burning during the CNO-flash. 



third dredge up episodes in low-mass AGB stars is of special 
importance for the scenario outlined in the previous para- 
graph, as third dredge up erodes the thin He-buffer formed 
just after a thermal pulse, ft is worth noting that mass-loss 
rates during the RGB and AGB phases have been adopted 
acc ording to the prescri p tions of ISchroder fc Cunta l|2005t ) 
and IVassiliadis fc Woodl (|l993h . This choice, together with 
the absence of third dredge up, is expected to produce a re- 
alistic initial final mass relation. The only exception is mass 
loss just before the final helium flash, which was tuned to 
obtain an AFTP/LTP episode according to the scenario pro- 
posed here. It is worth mentioning that standard AGB winds 
has been adopted during the He-burning phase after the last 
helium flash (AFTP/LTP), when the star becomes a born 
again AGB star. No mass loss has been adopted during the 
high temperature (T e ff ) phase of pre-white dwarf evolution. 



4 RESULTS 

In the previous section, we mentioned that several sequences 
have been calculated for the present study. As all our DIN se- 
quences display similar characteristics, we will first describe 
the evolution of a choosen reference sequence. Then, we will 
discuss how differences in progenitor/remnant masses, initial 
metallicities and He-buffer thicknesses may affect the course 
of evolution. In particular, we will show in later subsections 
the range in the predicted abundances and timescales of all 



our computed sequences, as well as the characteristics of 
those sequences that did not undergo DIN episodes. 



4.1 Detailed description of a reference sequence 

The evolution of our reference sequence in the HR-diagram is 
shown in Fig. [2] It corresponds to an initially M — 0.85 Mq 
(ZAMS) model star with Z = 0.0001. After evolving through 
the central H- and He- burning, the star starts climbing the 
early AGB. This is so until the development of the first 
thermal pulse (t agc = 11757.6 Myr, M = 0.584M©), which 
marks the beginning of the TP- AGB phase. The mass of the 
hydrogen-free core at that point is Mhfc = 0.515MQ. Due 
to the low mass of the star at the beginning of the TP- AGB, 
and owing to the relatively high winds (M ~ 10 _7 Mo/yr) 
the star suffers during the TP-AGB and the long inter- 
pulse period (P ~ 530000 yr), the star experiences only 
3 thermal pulses. Mass loss during the very end of the TP- 
AGB was tuned for the last thermal pulse to develop when 
the star was just leaving the AGB (an AGB final thermal 
pulse; AFTP). The mass at the moment of the AFTP is 
M = 0.54403Mq with a hydrogen-free core mass of Mhfc = 
0.53839Mq . After winds have reduced the mass of the star 
to 0.53946M Q , the star leaves the AGB (T cff > 10000K) and 
evolves towards the white dwarf stage with a H-envelope of 
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Figure 3. Evolution of the chemical profiles of our reference se- 
quence (0.53946Mq). Panels A, B, and C correspond to the in- 
stants denoted by the same letters in Fig. [2] Dashed lines in the 
middle panel indicate profiles just after the development of the 
convective zone driven by the H-flash. 
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Figure 4. Predicted lightcurves for the DINs of remnants coming 
from an initially M = 0.85M Q , Z = 0.0001 star. Different rem- 
nant masses were obtained by altering the wind intensity during 
the AGB final thermal pulse. Visual magnitudes were obtained 
from the predicted luminosity of the sequences and adopting the 
bolometric corrections derived by Flower (1996). 



M™ E ' WD = 1.07 X 10~ 3 M Q @ which corresponds to a total 
H-content of M™ E - WD = 6.6 x 1(T 4 M Q , and a He-buffer 
thickness of only M^* 13 = 1.37x 1O~ 4 M . Before enter- 
ing the white dwarf stage H is reignited, this produces the 
kink in the HR-diagram (see Fig. [2j just before the point 
of maximum T e g. As a consequence of this reignition, both 
the mass of the H-envelope and the total H-content of the 
star are diminished to 8.5 x 10~ 4 M Q and 4.2 x 1O _4 M , 
respectively. Conversely, the mass of the He-buffer zone is 
increased up to MauScr = 4.38 x lO~ 4 A/0. At this point 
the sequence begins its cooling as a white dwarf with a thin 
He-buffer see (Fig. 0. The chemical profile at that point is 
shown in Fig. [3] (panel A). From that point onwards, the 
star cools and contracts as a white dwarf. Simultaneously, 
chemical profiles are altered by element diffusion. In partic- 
ular, H from the envelope and C from the intershell region, 
just below the He-buffer (see Fig. [1]), diffuse through the 
He-buffer due to chemical diffusion triggered by the steep 
chemical profiles left by the previous evolution. Fig. [3] (panel 
B) shows the chemical profile 7.58 x 10 6 yr after the white 
dwarf started its cooling evolution. At this point the amount 
of C and H diffused into the He-buffer and, consequently the 
rate of energy generated by CNO burning starts to increase. 
From this point onwards CNO burning increases slowly un- 
til log Lh I Lq ~ 1.5 when the rate of energy generated by 
CNO burning leads to the development of convection (about 
23000 yr after CNO burning started to increase). This con- 
vective zone drags H from regions with higher H mass frac- 
tions and thus leads to higher energy generation rates. As a 
consequence the star enters a quicker phase and in less than 
200 yr Lh rises to log Lh / Lq ~ 7.5. This sudden energy 
injection causes a very quick expansion of the layers on top 
of the H-burning shell and the star is balloned back to the 
giant branch in a few years (see Fig. [2]| . In Fig. |4]we show the 
predicted visual lightcurve for the outburst episode — t = 
denotes the instant of maximum energy release by the CNO- 
flash, panel C of Fig. [3] see also Fig. [2] Visual magnitudes 
have been obtained from the predicted luminosity of the se- 
quences and adopting the bolometric corrections derived by 
Flower (1996) for normal H-rich abundances. As can be seen 
after suddenly rising from My ~ 9 to My ~ 2 the star in- 
creases its visual magnitude much slower than a classical 
nova, at a speed of ~ l m /yr. Specifically, the brightening 
speed just before reaching its maximum (at My ~ —5.2, 
~ 8 years after the H-flash) is dMy/dt ~ 1.6™ /yr. Just af- 
ter the maximum energy release by the H-flash, when the 
star moves to low temperatures and luminosities, the con- 
vective region driven by the H-flash merges with the (very 
thin) external convective region caused by low tempera- 
ture opacity. Then, the material from the very outer layers 
of the star is mixed with the material of the flash-driven 
convective region. Consequently, the surface abundances of 
the model are altered. The amount of He rises significantly 
(Xjj.e = 0.61 by mass fraction) while the H mass fraction is 
reduced (Ah = 0.39). Then, the DIN gives rise to an object 
with a mild H-deficiency. Also N is increased at the surface of 
the star. As the inner regions, which have experienced CNO 
burning, are N-enhanced, the surface abundance of N during 



2 We define the H-envelope as those regions with a H-content 
greater than 1 X 10 — 4 by mass fraction. 
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Figure 5. Evolution of the total H-content of our reference white 
dwarf sequence. The solid line shows the evolution in the case that 
AGB-like winds are considered during the second AGB phase af- 
ter the CNO-flash, while the dashed line indicates the evolution 
of the H-content in the case no AGB-like winds have been con- 
sidered. Note that regardless the intensity of winds in this stage 
the final H-content of the DA white dwarf is not modified and of 
the order of M H ~ 5 x 10~ 5 M Q 



the outburst increases to 2.2 x 10 -4 , significantly higher than 
the mass fractions of C and O (1 x 10~ J and 1.1 x 10~ 6 , re- 
spectively), and higher than the surface N abundance before 
the outburst (N ~ 1.56 x 10~ 5 ). Once back at low effective 
temperatures the star is expected to undergo strong winds 
similar to those of AGB stars. The time spent at low effective 
temperatures during this stage will be strongly dependent 
on mass loss assumptions. If standard AGB mass loss is as- 
sumed then the star reheates back to T e g J> 10 000K in about 
~ 50 yr, while in the case of no significant mass loss the star 
stays at low effective temperature for about ~ 260 yr. Un- 
der the assumption of standard AGB mass loss, the total 
mass lost by the object before it starts to contract again 
is of 1.6 x lO _4 M0(see Fig. [SJ. Finally the star returns to 
the white dwarf cooling track. The model returns to its pre- 
outburst stage (log L/Lq ~ —0.25 and log T c s ~ 4.6) in only 
~ 10 yr, much faster than in its pre-outburst cooling, which 
took ~ 7.6 x 10 6 yr. This should not be a surprise as only the 
outer layers are affected by the eruption and, so, the inte- 
rior of the white dwarf still resembles the structure prior to 
the eruption (see Fig. [6]). Only the outer layers, which have 
much shorter thermal timescales, cool and contract during 
this second white dwarf stage. The main consequence of this 
fast cooling is that diffusion has not enough time to alter the 
composition of the inner regions of our new white dwarf at 
high effective temperatures and the luminosity of the CNO- 
burning shell remains very low. Then, no new CNO-flash 
does take place this second time. 

4.2 Dependence on physical and numerical 
parameters 

In order to understand to what extent the outburst proper- 
ties depend on the remnant properties, we have performed 
several simulations under different assumptions. Specifically, 
we have calculated the outburst for remnants of different 



masses and also for remnants coming from progenitors of 
different masses. We also explore how thin the He-buffer 
needs to be in order to obtain a DIN. The latter is an im- 
portant issue to understand how frequent these events might 
be. 

To assess the importance of the He-buffer thickness for 
the development of DIN, we have computed four additional 
sequences by shutting down winds in our reference sequence 
at different times immediately after the maximum of He- 
burning, during the final thermal pulse. Then, the remnant 
leaves the AGB with thicker H-rich envelopes (and slightly 
different final masses). As a consequence, when these rem- 
nants contract towards the white dwarf stage, the reignition 
of H leads to thicker He-buffers. Hence, we are able to follow 
the cooling of white dwarfs with different He-buffer thick- 
nesses (see Table [I]). Interestingly enough, the H-envelope 
mass of the white dwarf is mostly dependent on the metal- 
licity and not on the detailed treatment of mass loss. 

As shown in Table Q] all the remnants derived from our 
0.85Mq Z = 0.0001 sequence, reach the white dwarf stage 
with very similar H-rich envelopes (within a ~ 10% differ- 
ence in mass). This is because the H-rich envelope of white 
dwarfs is mostly dependent on the mass of the remnant and 
initial metallicity. On the other hand, remnants display sig- 
nificantly different He-buffers thicknesses (with differences 
up to a factor 5). As a consequence of this difference, the 
CNO-flash developes at different points (luminosities) in the 
white dwarf cooling track. In fact, the thicker the He-buffer, 
the longer time needed for diffusion to mix enough H and 
C in the He- buffer and trigger the CNO-flash (see Table [2]). 
Also, we find that in our 0.85 Mq Z — 0.0001 sequence, if the 
He-buffer mass is AfBuffer > 2.1 x 1O _3 M0 the sequence does 
not experience a DIN. In Fig.[4]we show the lightcurves com- 
puted for these sequences. In Table Owe show the expected 
abundances for major chemical elements at the surface of the 
star during the nova outburst. As mentioned in the previous 
section, the mixing of material from the He-buffer with that 
at the H-rich envelope by convection, leads to an enhance- 
ment of the surface abundances of He and N at the expenses 
of the H abundance. Depending on the relative masses of the 
He-buffer and the H-rich envelope the outcome of such mix- 
ing gives different values of H, He and N. In fact, for those 
sequences with thicker He-buffers, the amount of He mixed 
during the CNO-flash is higher. This can be clearly appreci- 
ated from both Tables [T] and [2] the remnants that enter the 
white dwarf stage with relatively thicker He-buffers are those 
who end with larger He-abundances after the eruption. In all 
the cases, however, only a mild (log Nh/Nhc ~ —0.15. ..0.6) 
He-enrichment is attained. It is worth mentioning that the 
amount of H burnt during the nova event is not significant 
and the star will ultimately end as a H-rich white dwarf (DA 
white dwarf). 

In order to explore to which extent the mass of the 
remnant and its previous evolution would be relevant for 
the properties of DIN events, we have calculated additional 
sequences for different progenitor stellar masses. As in the 
case of our reference sequence, we have evolved these se- 
quences from the ZAMS up to the AGB phase and tuned 
mass loss during the last thermal pulse in order to obtain 
a final thermal pulse when the H-envelope mass was al- 
ready very thin (i.e. an AFTP/LTP). Again, as it happens 
with our reference sequence, (standard AGB) winds during 
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unchanged. 



the AFTP/LTP erode a significant fraction of the remain- 
ing H-envelope, inhibiting an important reignition of the H- 
burning layer. Thus, white dwarfs with very thin He-buffers 
are obtained. Namely, we have performed 3 additional cra- 
dle to grave sequences with Z = 0.0001 and ZAMS masses 
of M = 1M©; 1.25M and 1.8M . These sequences end as 
white dwarfs with thin-He buffers and masses of 0.55156Mq 
(M h = 2.1xlO- 4 M Q ; M Buffcr = 3.6xlO~ 3 M ), 0.59606M Q 
(M H = 2.8 x 1O _4 M ; M Buffcr = 9.8 x 10 _5 M Q ) and 
O.7O492M (M H = 1.8 x 1O _5 M ; M Bu ff cr = 5.9 x 1O _5 M ). 
While the first two sequences suffered a DIN (see Table[5]for 
their outburst properties), the last (more massive) one did 
not. In fact in order to explore how robust this result was, we 
have followed a similar approach than we did with our ref- 
erence sequence. By tuning the amount of mass lost during 
the very last AGB phase, we obtained white dwarfs with 
very similar masses but significantly different He-buffers 
and H-contents, namely: O.7O492M (M H = 1.8 x 1O~ 5 M ; 
M Bu ffcr = 5.9xlO~ 5 M ), O.7O496M (M H = 4.3xlO- 5 M ; 
M BuSel = 5.9xlO _5 M ), O.7O5OOM (M H = 7.5xlO- 5 M ; 
M BuSel = 8 x 1O _5 M ), O.7O54OM (M H = 6.2 x 1O- 5 M ; 
MbuBct = 4.8 x 1O _4 M ). In any of these sequences did 
we obtain a DIN. Thus it seems that DINs do not oc- 



cur in white dwarfs with masses significantly higher than 
the canonical O.6M value. To explore further the impor- 
tance of the white dwarf mass for the development of DINs 
we computed additional sequences by artificially stripped 
mass during the fifth thermal pulse of our 1.8M sequenc^f] 
to obtain a O.62361M white dwarf with a thin He-buffer 
(M H = 2.6 x 1O _4 M ; M Bu ff cr = 7.5 x 1O _5 M ). This se- 
quence did experience a DIN and its properties are described 
in Table H 

For the sake of completeness we have also performed 
two additional full sequences with progenitors of M = 1M 
but two different metallicities; Z — 0.001 and Z=0.01. These 
sequences end as white dwarfs of O.558O9M (Mh = 2.9 x 
1O- 4 M ; M Buf f er = 2.8 x 1O _4 M ) and O.52382M (M H = 
2.5 x 1O" 4 M ; M Bu ffor = 5.5 x 10~ 4 A/ Q ). The first sequence 
did experience a nova event (see Table [2]), while the second 
one did not. 

Although the exploration of the parameter space we 
have made is not completely exhaustive, we can draw some 
preliminary conclusions. To begin with, it seems that the 

3 Under standard AGB winds the sequence experiences 17 ther- 
mal pulses before ending as a O.7O492Af white dwarf. 
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occurrence of DINs is confined to low-mass, low-metallicity 
remnants. This is true regardless the fact that our most 
massive remnant, 0.70492Mq , had a very thin He-buffer 
(M B uff or = 5.9 x 10~ 5 M©). The absense of DINs in higher 
mass or metallicity remnants with thin He-buffers is cer- 
tainly not only because of the lower H-content that higher 
mass/metallicity white dwarfs are supposed to have. In fact, 
the O.52382M0 Z=0.0I sequence does not show significant 
differences neither in mass nor in the values of Mu and 
M BuBcv from that of our reference 0.53946 M© Z=0.0001 se- 
quence. The key point that seems to prevent the DIN event 
in our more metallic sequence is related to the luminosity 
of the CO-core. In fact, despite the slightly lower mass and 
slightly higher carbon content, our O.52382M0 Z=0.0I se- 
quence cools significantly faster (« 50% at logT c ff = 4.6) 
than our reference 0.53946 Mq Z=0.000I sequence. Con- 
sequently, the core luminosity is roughly ~ 50% higher 
than that of the reference sequence. This, coupled to the 
fact that the CNO-burning shell is « 25% less luminous 
(at log T^ff = 4.6) than the reference sequence case, makes 
CNO-burning stable (see section[4]3]) and thus no DIN devel- 
opes. The faster evol ution of higher metal licity warm white 
dwarfs was shown in iRenedo et all (|2010l ). It seems, then, 
that higher metallicity tends to prevent CNO-flashes, mostly 
due to a faster evolution during the warm white dwarf stage 
which leads to a somewhat higher gravothermal energy re- 
lease in the core. In the case of our highest mass white dwarf 
(~ 0.70500M Q ), the absense of CNO flashes is related to the 
low luminosity of the CNO-burning shell, which is, at most 
(at logTefi ~ 4.55), Lcno ~ L COTE (see section |4~3|) . Then, 
our numerical experiments, suggest that higher metallicity 
or masses tend to prevent the occurence of DIN events due 
to the stabilizing effect of a low Lcno /L cole value (see sec- 
tion [OJ. 

It must be noted that the final masses of our white dwarf 
remnants are in agreement with our present understand- 
ing of the ini tial-final-mass rela tionship for low/intermediate 
mass stars (|Salaris et al.l l2009h . and thus we expect our 
white dwarf models to be physically sounding. In fig 
we show the lightcurves and effective temperature evolution 
during the outburst for sequences with different progenitors, 
metallicities and final white dwarf masses. 

Finally, our simulations show that DINs have charac- 
teristic properties that make them easy to recognize. First 
their visual lightcurve displays a sudden increase, in less 
than a month, of 6 to 7 magnitudes followed by slower in- 
crease in which it increases from My ~ 2 m to its maximum 
of My ~ — 5.5 m in a few years — ~ 2 and ~ 10 yr depend- 
ing on the remnant properties. Its brightening speed just 
before reaching their maximum would be of I to 4 magni- 
tudes per year. Also, as the star brightness increases the 
star should become cooler and yellowish (see Fig. 0. Such 
lightcurve w ould be probably classified as an extremely slow 
rising nova (|lben fc MacDonaldl Il986h . Secondly, the out- 
burst abundances (see Table[3]) of all our sequences are char- 
acterized by log Nh/Nh c ~ —0.15. ..0.6 and high N abun- 
dances of 10" 4 - 10" 3 and C>0, by mass fraction. 

4.3 The stability of CNO-burning on white dwarfs 

In this section, we will try to clarify the reasons beneath the 
instability of the CNO-burning shell by looking for a quan- 
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Figure 7. Predicted lightcurves for the DINs of remnants with 
different masses and descending from different progenitors. See 
text for details. Visual magnitudes were obtained from the pre- 
dicted luminosity of t he sequences and adopting the bolometric 
corrections derived by iFlowe 
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Figure 8. Evolution of the value of the stability parameter a 
for selected white dwarf sequences. Note the significantly lower 
a- values of the high mass (O.7O5OOM0) and high metallicity 
(O.52382Af , Z=0.01) sequences which denote the stability of 
CNO-burning for those sequences. Also note the very low cr-value 
of our reference sequence during its second evolution of the white 
dwarf stage. 



titative measure of the stability of H-shell burning. In this 
co nnection, we have closely followed the criterion presented 
bv lYoon et alj (|2004T ) and extended it to be used with mild 
bur ning shells in white dwarfs. 

lYoon et al.l (|2004l ) analysed the stability of a burning 
shell by investigating the reaction of its mean properties 
under a uniform expansion of the shell and the surrounding 
envelope. Specifically, they defined the burning shell (where 
mean quantities are calculated) as those layers where the en- 
ergy generation rate is e C NO > ecno at shell pca k/0.002. By 
assuming that the inner boundary of the burning shell (ro) 
remains constant, neglecting the gravitational energy release 
of the unperturbed model, and neglecting both the lumi- 
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nosity of the core (lo = 0) and its perturbation (Slo = 0), 
lYoon et alj |2004 ) found that the temperature perturbation 
9 = 5T/T follows the equation 



othin ~ £t x + K T y — 4 y + a T (e p x + k p y) . 



(4) 



r th e = a8. 



(1) 



Where Tth is the thermal time scale of the shell and a is 
given by 



er + k t - 4 + (e P + a s + K p ) 



(2) 



Where er, e p , kj, k p are the logarithmic derivatives of 
e(T, p) and «(X, p) while ai, ap are the logarithmic deriva- 
tives of p(T, P). a s is a parameter describing the geom etrical 
thickness of the burning shell (see lYoon et al.1 12004! ). Then, 
a is the parameter quantifying the shell instability, so that 
a < (a > 0) corresponds to stable (unstable) burning. 

In our case, an important shortcoming of this approach 
comes from the fact that the CNO-burning shell is not 
clearly detached from the outer pp-burning shell. Thus, the 
choice of the borders of the simplifie d burning shell i s not 
as easy as in usual burning shells (see lYoon et aUl2004 ). For 
the present analysis we have choosen the lower boundary of 
the shell as the point at which e C NO = ccno at shell pca k/10, 
which is close to the point at which the structural effects 
of the burning shell start to be apparent. For the outer 
boundary we have choosen it to be the point closest to 
the shell maximum at which either ecNO — e pp or ccno — 
£CNO at shell pcak/10. The other importan t issue comes from 
the fact that the core luminosity lo (see lYoon et al.l |2004 ) 
can not be neglected in our case, due to the relative im- 
portance of the gravothermal energy release from the core 
as compared to the shell luminosity (?cno = Is — lo, where 
i s stands for the luminosity at the outer boundary of the 
burning shell). Also, the perturbation of the luminosity at 
the inner boundary of the burning shell (<5Zo) ca nnot be ne- 
;lecte d eitheiQ- Within the scheme presented by lYoon et al.l 
20041 ) Sl can not be derived and, thus, we will consider it 
to be proportional to Sl 3 . It is clear that, for the reasons 
mentioned above, we can not expect from this approach to 
get a quantitative prediction of the instability of the CNO- 
burning shell, but we will show that its qualitative descrip- 
tion is very useful to understand the causes beneath the shell 
instability. 

Under the assumptions and definitio ns mentioned above 
it is easy to show that the criterion of lYoon et al.l (|2004 
is modified so that the parameter a quantifying the shell 
instability is defined as 



e T x + ti T y - 4j/ + 



(ctpa s — 1) 



jr (e p x + a a y + k p y) 



'ad 



(3) 



(otp ot 3 —1) 



The main differ ence with the criterion presented by 
lYoon et al.1 (|2004 ) comes from the presence of the param- 
eters x — (l a — lo)/l a and y — (8l 3 — Slo)/5l a . For the 
present analysis we h ave adopted extre me values for y; 
y = 1 (5l = 0, as in lYoon et al.1 (12004 )) and y = 2 (i.e. 
Slo = —Sis). 

In all our sequences ot s < 0.16 so the behaviour of the 
burning shell is not far from that of thin burning shells, 
where 

4 Note that this is also true for those cases in which Iq ~ 0. 



In the absence of a core luminosity (x = 1) 0thin> and thus 
instability, is dominated by er - However, in our sequences we 
find the factor x to play a very important role as its values 
range from x ~ 0.5 for the 0.53946 MQ-sequence just before 
the flash to less than 0.2 for the 0.70500Mq -sequence at its 
maximum value. An exploration of our sequences shows that 
the term er x is dominant in the behaviour of a. 

In Fig. [8] we show the evolution of the parameter a for 
some selected sequences. As expected, the quantitative pre- 
diction of the instability is not reproduced. However quali- 
tatively, Fig. [8] shows that instabilities develope during the 
rise in the parameter a which is caused by the increase in er 
due to the cooling of the shell. However, this increase in er is 
counterbalanced by the decrease in a; as a consequence of the 
dimming of the shell. The interplay between this two main 
factors is crucial for the development of the instability. Dif- 
fusion leads to an increase of the CNO-burning shell power, 
which increases x and might trigger the flash. If, by the time 
diffusion starts to change the H-content at the H-burning 
layers the shell is already very dim H-burning remains stable 
and no DIN develops. The existence of a thicker He-buffer 
delays the effect of diffusion and thus, if the buffer is thick 
enough H-burning remains stable during the white dwarf 
cooling. From Fig.|5]we can also understand why higher mass 
sequences (M ~ 0.70500M Q ) do not undergo DIN events de- 
spite the very thin helium buffers. The reason is that burn- 
ing shells are too dim as compared with the core luminosity 
and, thus, very low x values render H-burning stable. In the 
case of the high metallicity (Z=0.01) 0.52382Mq -sequence 
discussed in previous sections, its faster cooling (than our 
reference Z=0.0001 sequences of similar mass) is tied to a 
higher core luminosity which makes CNO-burning stable. 

Thus, our analysis suggests that CNO-flashes depend 
on the intensity of the CNO-burning shell (compared to the 
core luminosity, i.e. x) and its temperature (i.e. ex)- Both 
factors depend on diffusion velocities, which mixes H into 
deeper, hotter, layers increasing the CNO-burning power, 
and the thickness of the He-buffer which delays the contact 
between C and H. Finally, the relative intensity of the CNO- 
burning shell compared to the core luminosity depends on 
the mass and (progenitor) metallicity of the white dwarf, a 
fact that seems to prevent DIN events in higher mass and 
metallicity remnants. 



4.4 Effect of extramixing in the CNO-flash driven 
convective zone 

In the interests of completeness we have also analysed the 
impact of extramixing in the flash-driven convective zone. 
In this connection, we have recalculated the DIN event for 
our reference sequence including extramixing beyond con- 
vective boundaries durin g the flash. Th i s was done within 
the picture developed by iHerwig et al.l (|l997l ) of an expo- 
nentially decaying overshooting (OV) with / = 0.016 (see 
IHerwig et al.lll997l for a definition). It is worth noting that, 
although the physical picture behind this prescription might 
not be correct, it can be considered as an effective way to 
include extramixing processes that take place at conv ective 
boundaries (e.g. gravity waves; see IHerwig et alj|2007l ). 

The main consequence of the inclusion of OV during 
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the development of the CNO-flash convective zone is that H 
is dragged deeper, into hotter regions of the star. This en- 
hances the CNO-burning luminosity during the flash, which 
can be as high as Lcno ~ 1O 1O -L , significantly higher than 
the value of Z/cno ~ 1O 7 L predicted by no-OV sequences. 
Also, due to the fact that the convective zone reaches deeper 
during the flash, much more material from the N-rich He- 
buffer (and even from the C-rich intershell) is dragged to the 
surface during the outburst. As a consequence the surface 
abundances during the outburst are significantly enhanced 
in He and N (see Table[3]). Also, surface C and O are strongly 
increased to ~ 10 -3 (by mass fractions), as compared with 
the 10~ 5 — 1CP 6 resulting from the case when OV is ne- 
glected. It is worth mentioning that the increase in C and O 
would be higher if the intershell abundance of the model is 
richer in C and O (see section [5T2l and Table |3). Also, as the 
CNO energy release is orders of magnitude higher when OV 
is considered, the event is more violent, reaching the giant 
stage in less than a year and showing a much steeper out- 
burst lightcurve (see Table[2]). Also, when OV is included the 
duration of the giant stage (ts, Table 0) is also enhanced, 
being of the order of ~ 10 4 yr. This is because much more 
H is consumed in this case and consequently much more en- 
ergy is released, allowing for a longer giant stage. Also as a 
consequence of this higher H-burning, the amount of H re- 
maining in the star by the time it reenters the white dwarf 
stage is about Mh ~ 1O _5 M , six times smaller than in the 
no-OV case (M H ~ 5.7 x 1O~ 5 M J| 

4.5 Comparison with llben fc MacDonaldlll986l 

In this sectio n we compare our results w ith those of the DIN 
calculated bv llben fc MacDonaldl (|1986). In Table[2]we show 
the main characteris tics of our simulations tog hether with 
those ext racted from Iben fc M acDonaldl (|l986h . An exami- 
nation of llben fc MacDonaldl (|l986l ) results show that all the 
predicted quantities fall within the range predicted by our 
simulations. In particular, it is worth noting that the time, 
luminosity and temperature at which the nova takes place in 
the simulation of Iben & MacDonald (1986; IM86) are very 
similar to those of our sequence of similar metallicity (al- 
though different mass). The same is true for the abundances 
of He and H during the outburst. In vi ew of the many nu- 
meric al and physical differences between llben fc MacDonaldl 
(1986) and our simulations, we consider the similarities of 
both results as a strong indication for the robustness of the 
results and the feasibility of the scenario. In particular, we 
feel confident that the results presented in Tables [2] and [3] 

reflect the main characteris tics of DIN outburts . 

Regarding the work bv llben fc MacDonaldl (|l986l ) some 
comments need to be done. The authors conclude that if 
the He-buffer thickness is smaller than a critical value given 
by M™/er ~ AMh/10, where AM H stands for the in- 
crease of the H-free core during the interpulse phase on 
the AGB, then a DIN occurs. This c o nclus ion has been 
recently used by lLawlor fc MacDonaldl l|2006r i to conclude 
that DINs can not occur in white dwarfs descending from 



progenitors of masses lower than 3M . This result is com- 
pletely at variance with the results presented in the previ- 
ous sections where all DINs come from low mass progen- 
itors. However, it should be mentioned that the condition 
^Buffer ~ AMh/W seems to be based on a single se- 
quence, and thus it may be wrong to extrapolate it to rem- 
nants of different masses. During the present study we have 
not found such "critical value" to be fulfilled. In fact, in 
our 0.85 M sequences we find the frontier between the de- 
velopment or not of a diffusion nova to take place some- 
where between our 0.54148M Q (M Buffer = 2.3 x 1O _3 M ) 
and our 0.54148M© (M Bu ffer = 2.06 x 1O" 3 M ) remnant. 
For these sequences, the increase of the H-free core be- 
tw een thermal pulses i s AM h = O.OIOIMq, very similar 
to llben fc MacDonaldl |l986) value. Then, the critical He- 
buffer thickness is M B ^jj er ~2x 10 ~ 3 ~ AM H /5. This 
is two times larger than suggested by llben fc MacDonaldl 
(1986). On the other hand, for our 1.8M progenitor we 
find that while our O.6236M (M Bu ffcr = 7.5 x 1O~ 5 M ) 
white dwarf remnant does experience a CNO-flash, our 
slightly more massive remnant of O.6239M (MBuffer = 
3 x 1O _4 M ) does not. Then, as the increase of the H-free 
Core for the 1.8M sequence during the fifth thermal pulse 
is AMh = 0.009, then the critical He-buffer thickness is 
M Bu}jer ~ AM H /a with a a number between 30 and 120. 
Finally, let us mention that when we followed the AGB evo- 
lution of our 1.8M sequence up to 17 thermal pulse, the 
growth of the H-free core during the interpulse phases is 
AMh = O.OO59M . For the white dwarf remnants created 
from this progenitor, even in the case with the thinner He- 
buffer, 0.70492M© (MBuffer = 5.9 x 1O" 5 M ), we did not ob- 
tain a DIN. Even though in this case MBuff cr = AMh/100, 
no DIN developed. This is in line with our analysis of Sec- 
tion l4.3l that suggests that for these white dwarf masses DIN 
events will not develope, regardless of the He-buffer mass. 
Then, we conclude that the previously mentioned critical 
He-buffer mass for the development of a DIN is misleading. 
As mentioned in section 14.31 it is clear that the He-buffer 
thickness is not the only ingredient that determines whether 

a DIN will oc cur. 

Finally, llben fc MacDonaldl (|l986l , speculattfl that 
DINs might be recurrent. As mentioned in the introduc- 
tion, th e occurence of recurrent DIN episodes has been in- 
voqued (|D'Antona fc Mazzitellilll990T ) as a possible mecha- 
nism to explain the existence of DA white dwarfs with thin 
H-envelopes (as those inferred by some asteroseismological 
studies; e.g. ICastanheira fc Keplerll2009] ). Our computation 
of the aftermath of the DIN, shows that the star evolves 
much faster during its second cooling stage as a white dwarf 
and thus, diffusion has no time to significantly alter the 
chemical profiles of the inner regions of the white dwarf. 
As a consequence, H-burning during the second white dwarf 
evolution remains stable and no recurrent DIN events can 
develope. 



5 Note that this factor corresponds also to the difference between 
the duration of the giant stage in the no-OV case, T3 = 260 yr, 
and the OV case T3 = 1367 yr. 



° Due to numerical difficulties they stopped their simulation once 
the star returned back to a giant configuration. 
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5 POSSIBLE OBSERVATIONAL 
COUNTERPARTS 

5.1 Speculations on the origin of CK Vul 

In the previous sections we have presented the predictions 
of the DIN scenario. As these events have been suggested to 
look like "very slow novas" it seems natural to compare their 
predictions with some of those objects. The most famous 
and enigmatic one is Hevelius nova CK Vul l|Shara et al.l 
1 19851 ). Here, we will compare the predictions of the DIN with 
observations from CK Vul and with the predictions of the 
best available scenario, the very late thermal pulse (VLT P; 
lHarrisorJ[l996| , lEvans et al.ll2002l and lHaiduk et alj|2007h . 

5.1.0.1 Surface composition, H-abundance One of 

the most remarkable difference between VLTP models and 
observations of CK Vul come from the significant H emis- 
sion detected in the ejected material of CK Vul, in contrast 
with the extreme H-deficiency detected in the ejected ma- 
terial of t he two bonafide VL TP objects; V4334 Sgr and 
V605 Aql (|Haiduk et al.ll2007T ). In this connection, the DIN 
scenario predicts a high H-content in the envelope during 
the ourburst and, thus of the ejected material (see Table 
H), contrary to the predictions o f the VLTP scenario (see 
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that the most abundant CNO element predicted by the DIN 
scenario is N, which is clearly det ected in CK Vul surround- 
ing material (|Haiduk et alj|2007h . 

5.1.0.2 Absence of an old planetary n ebula From 
the lo w dust emission surrounding CK Vul, lHaiduk et al.l 
conclude that CK Vul must be an evolved object 

10 5 — 10 6 yr older than the last time it experienced steady 
high mass loss (e.g. on the TP-AGB). This time is orders of 
magnitude longer than the time from the departure of the 
TP-AGB to a VLTP (~ 10 4 yr). In fact both V4334 Sgr and 
V605 Aql show old planetary nebulae sur rounding the new 
ejected material |Kimeswenger et al.l2007T ). On the contrary, 
the DIN scenario predicts a much longer timescale between 
the cesation of high mass loss and the time of the CNO- 
flash. As shown in Table [2] this timescale is of the order of 

10 6 — 10 7 yr, in closer agreement with the observations of 
CK Vul. 

5.1.0.3 Ejected Mass Due to the fact that the expand- 
ing layers in a DIN episode are those above the region of the 
CNO-flash, DIN models have a very stringent upper limit to 
the mass that can be expelled during the outburst. As an ex- 
ample, the location of the H-flash in our reference sequence 
is < 1O _3 M0 below the surface and the expanded envelope 
once the star is back on the AGB is of about ~ 1O _4 M0. 
Thus, we can safetly conclude that the ejected mass in the 
DIN scenario will be of ~ 1O _4 M0 — and surely smaller 
than ~ 10 _3 Mm. These values are within the inferred masses 



for the ejected material in CK Vul, which is of 10 M ( 



s 



5 x 1 2 Mq depending on the assumptions (|Haiduk et al.l 
l2007h . 

5.1.0.4 Outburst lightcurve CK Vul was discovered 
at its maximum brightness on the 20th of June of 1670 as 
3rd magnitude star. Although no description of the out- 



burst lightcurve exist, IShara et al.l (| 19851 ') conclude that it 
is very unlikely that the star was visible in the previous 
season. Thus, only a lower limit to the brightness speed of 
~ 2 magnitudes per year can be established. This limit is 
within the values predicted by the DIN scenario (see Ta- 
ble [5J). After reaching maximum brightness it faded and re- 
covered twice before finally disappearing from view. These 
fadings are remarkably si milar to those observe d in born 
agai n AGB stars (V605 Aq l; lDuerbeck et alfcool and V433 4 
Sgr; iDuerbeck et alj|200oh and RCrB stars l |Clavtor]|l996h . 
In fact, this particular lightcurve was one of the main ar- 
guments to link CK V ul to the VLTP scenario l|Harrisonl 
1 19961 lEvans et alJl2002T ). RCrB-like fading are usually asso- 
ciated with carbon dust so one should wonder if such fading 
episodes could take place in atmospheres which are not C- 
rich, as those predicted by our DIN simulations (see Table 
[2|. In this connection, we find suggestive tha t there is one 
"Hot RCrB", MV Sgr |De Marco et all 120021 1. which has a 
carbon and nitrogen composition very similar (C ~ 5 x 10~ 4 
and N ~ 9 x 10 -4 ) to those predicted by the DIN scenario 
(N>00 and N ~ 10~ 3 - 10~ 4 , see Table©, although MV 
Sgr shows a much stronger He-enrichment. We consider that 
MV Sgr suggests that DIN events might also display the fad- 
ings observed in RCrB like lightcurves, although we consider 
this just an speculation and should be analysed in detail. 
From these considerations we conclude that our present un- 
derstanding of the DIN-scenario is not in contradiction with 
CK Vul eruption lightcurve. 

5.1.0.5 Maximum brightness The maximum bright- 
ness is the main discrepancy between observations and mod- 
els (both for DIN or VLTP models). While observations of 

CK Vul suggest at the preferred distance M v CKVul ' 8, 

DIN models predict significantly lower values of Mv Um ~ 
—5.5. It is worth noting, however that this is even worst 
for VLTP models, in which fast- VLTP sequences predict 
My VLTP ~ —4 (according to the predicted luminosities in 
iMiller Bertolami fc Althausll2007h . 

5.1.0.6 Actual luminosity of CK Vul One of the 

main discrepancies between inferences from CK Vul and 
both DIN and VLTP scenarios comes from the radio 
emission. In fact, from the radio emission detected by 
lHaiduk et ail l|2007h these authors derive a very low lumi- 
nosity of ~ Lq for the ionizing object. This is done by con- 
sidering that hydrogen is ionized by the photons emited by 
the central object above the Lyman limit. Neither the DIN 
nor VLTP events can account such very low luminosities 
only 330 yr after the outburst. However, it is worth noting 
that such very low luminosity is derived from the assump- 
tion that the temperature of the central object is close to the 
effective temperature at which the number of ionizing pho- 
tons is maximum. From the data presente d by Haiduk et al.l 
l|2007h and adopting the expression from IWilcotsI l|l994h . 

F c = 1.761 x 10 48 x a(v,T e ) x v 0A x T e ~ 0AS x S x D 2 , (5) 

we derive that the number of ionizing photons is iV 7 ~ 
4 x 10 43 s _1 . As the number of ionizing photons is steeply 
dependent on the effective temperature of the star, such 
number of photons can also be attained by a star of T e g ~ 
10000A and L ~ IOOOOLq. These values put the central ob- 
ject back in the AGB region of the HR-diagram. Hence we 
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are tempted to speculate that the radio emission could be 
compatible with an object reheating back from the second 
stage in the AGB and thus, compatible with both VLTP 
and DIN scenarios. However, increasing the luminosity of 
the st ar should increase the IR-flux from the obscuring ma- 
terial ijEvans et all [2002). One may thus wonder if the IR- 
flux would be much higher than the observational limits (e.g. 
S12 ~mJv; lHaiduk et al.ll2007h . Some rough numbers show 
that this would be strongly dependent on the geometrical 
and physical assumptions on the obscuring material. For ex- 
ample, if we assume the star being obscured by a circumstel- 
lar shell located at r ~ 30AU (and assumed to behave as a 
blackbody), we obtain a shell temperature of Td us t ~ 720K 
and a flux of S12 ~ 1000 Jy, several orders of magnitud e 
above the observational upper limit (|Haiduk et al.l l2007f ) . 
On the contrary if the star is obscured by the presence of 
a pre-existing edge-on disk (n nner = 1AU, r out cr = 60 AU 
and height h — 1AU), as suggested by the ejection almost 
perpendicular to the line of sight, then the IR flux would 
be much lower. In deed, by assu ming a classical temperature 
profile (T oc r~ 3/4 ; lNattalll993l ') and a black body behaviour 
we (roughly) estimate S1 2 ~ 25 mjy, much closer to the 
observational constraints (|Haiduk et al.|[2007l ). 

We conclude that, if CK Vul is actually obscured by 
a thick disk, the radio emission could be explained by a 
cool reborn giant star (T off ~ 100007f and L ~ 10000L Q ) 
without an strong IR emission. This would help to reconcile 
the predictions of both VLTP and DIN scenarios with the 
observed properties of CK Vul. Then, a detailed modeling 
of the dust surrouding CK Vul seems to be key in order to 
reject a possible luminous cool central star as the present 
state for CK Vul. 

5.2 Possible connection with PB8 

As mentioned in Section 14.41 the inclusion of extramixing 
processes at the boundaries of the CNO-flash driven con- 
vective zone leads to outburst abundances rich in He, H 
and N (see Table Also, the inclusion of extramixing 
processes leads to a more violent flash and to an increase 
in the time spent as a giant (T3; see Table [2]) after the 
DIN event. During this short giant stage the reactivation 
of AGB winds is expected. Depending on the amount of 
mass lost during this stage, it could be possible that the 
new ejected material might form a new young planetary 
nebula once the star begins to contract and heat towards 
the second white dwarf phase. Such young planetary nebula 
would show very low dynamical ages and, together with its 
central star, a very atypical He, H and N rich composition. 
PB8 has been r ecently shown to have (qualitatively) simila r 
characteristics (|Todt et al.ll201ol . iGarcfa-Roias et al.| [2009'). 
namely low dynamical age of a few thousand years and a He- 
, H- and N-rich composition. Specifically, the central star of 
PB8 has been proposed as the prototype of [WN/WC], a 
new class of CSPN. In particular, the surface abunda nces of 
PB8-CSPN have been derived, bv lTodt et all (|2010H . to be 
[H/He/C/N/O]=[40/55/1.3/2/1.3]. These abundances can 
be accounted by our DIN sequences that include extramix- 
ing processes at convective boundaries — see Table[3] Indeed, 
as mentioned in Section [4.41 when overshooting is included 
at the CNO-flash convective zone, convection drags material 
from deeper layers of the white dwarf, reaching the upper 



layers of the C-rich intershell. Thus, the C and O surface 
abundances after a DIN event will be strongly dependent on 
the intershell composition of the star when departing from 
the AGB. In fact, if intershell abundanc es are assumed to 
be si milar to those of PG1159 stars (see IWerner fc Herwid 
2006) much higher C and O abundances are obtained af- 
ter a recalculation of the DIN event (see Table [3]). In fact, 
the abundances obtained in this case are in good q ualitative 
agreement with those derived bv lTodt et all (|2010l ) for PB8- 
CSPN ( i.e. He>H>N>C~ Q). Also, the dynamical age de- 
rived bv lTodt et all l |20ltf) . t = 2600 yr, is qualitatively sim- 
ilar to the giant phase duration of this sequence (T3 ~ 1600 
yr). Taking into account that our simulations correspond to 
a single, arbitrarily choosen sequence, we find the similarities 
with the observed properties of PB8 to be remarkable. Fur- 
ther exploration of this possible explanation for PB8 seems 
natural. 

6 FINAL REMARKS 

We have performed a detailed exploration of the diffusion- 
induced nova scenario (DIN) originally proposed by 
llben k, MacDonaldl (1 19861 ). Our calculations represent a sig- 
nificant improvement in our understanding of CNO-flashes 
in warm white dwarfs. 

Our main results can be summarized as follows: 

• We have identified a definite scenario leading to the 
formation of DA white dwarfs with thin He-buffers. Such 
white dwarfs are naturally formed in low-mass stars, that 
do not experience third dredge up during the TP- AGB, and 
suffer from either an AFTP or a LTP. 

• We have explored the parameter space of the DIN sce- 
nario and shown that there is a range of values of M*, 
Zzams and He-buffer masses for which DIN occur in phys- 
ically sounding white dwarf models. Our results suggest 
that DINs take place in white dwarfs with M* <; 0.6 and 
Zzams < 0.001 and thin He-buffers — as those provided by 
the scenario described above. 

• Our simulations provide a very detailed description of 
the events before, during and after the DIN event. In partic- 
ular, our results show that DIN events are not recurrent as 
previously speculated. Thus, DINs do not form H-deficient 
white dwarfs, nor DA white dwarfs with thin H-envelopes. 

• We have qualitatively described the mechanism by 
which the CNO-shell becomes unstable. Our analysis shows 
that the occurrence of CNO-flashes depends strongly on the 
intensity of the CNO-burning shell (as compared to the core 
luminosity), and its temperature. This seems to be in agree- 
ment with the fact that only our sequences with M* <J 0.6 
and Zzams < 0.001 and thin He-buffers experienced DIN 
events. 

• Regarding t he criterion presented by 
llben fe MacDonaldl (|l986r ) for the occurrence of DIN 
events, we find that such criterion is misleading as the 
He-buffer mass is not the only parameter that determines 
whether a DIN event will take place or not. In particular, 
for more massive remnants our simulations do not predict 
DIN for any possible He-buffer masses. 

• Our simulations provide a very detailed description of 
the expected surface abundances and lightcurves during the 
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outburts. In particular we find that typical lightcurves dis- 
play a maximum of My ~ —5.5, a brightness speed of 
a few magnitudes per year, and a mild He- and N- en- 
richment; with N ~ 10~ 4 — 10 -3 by mass fraction and 
log Nu/Nuc ~ —0.15. ..0.6. Also, in all our sequences we find 
surfaces abundaces with N > C > O by mass fraction. 

• We find that the inclusion of extramixing events at the 
boundaries of the CNO-flash driven convective zone leads to 
higher He, N, C and O abundances than in the case in which 
no extramixing is considered. Relative surface CNO abun- 
dances in these cases are N > C J> O (by mass fractions), 
although the precise values will be strongly dependent on 
the C and O composition of the He-, C- and O- rich inter- 
shell. 

• Finally, with the aid of our numerical simulations, we 
have discussed the possibility that Hevelius nova CK Vul 
and PB8-CSPN could be observational counterparts of DIN 
events. 

For CK Vul we find that despite discrepancies with ob- 
servations the DIN scenario offers a possible explanation for 
CK Vul which is as good as the best available explanation 
to date (the VLTP scenario). In particular the DIN sce- 
nario can easily explain the absense of and old planetary 
nebulae surrounding CK Vul and the presence of H in the 
ejected material, which cannot be understood easily within 
the VLTP scenario. Moreo ver, we suggest that the radio and 
IR emission from CK Vul (|Haiduk et alJl2007T ) could be ex- 
plained if the star is assumed to be in its giant stage and 
surrounded by a pre-existing dust disk. Such reinterpreta- 
tion of the observations would reconcile them with both the 
VLTP and DIN scenarios. A further exploration of this pos- 
sibility seems desirable. 

Regarding PB8, our simul ations show that the atypical 
abundances observed in PB8 dTodt et alJl201Ch and the low 
dynamical age of its nebula (|Garci'a-Roias et alj|2009h could 
be qualitatively understood within the DIN scenario if ex- 
tramixing is allowed at the boundaries of the CNO-flash 
convective zone. 
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Table 1. Properties of post-AGB remnants of an initially O.85M0 star with Z = 1 X 10 -4 which experiences a final thermal pulse (see text for details). M env is the mass of the outer 
layers with Xjj > 10 — 4 . A?Buffer Thickness of the almost pure He layers (He-buffer). Mjj total hydrogen content of the remnant. PRE— WD m( jicates values just after departing from 
the AGB, while the last three columns correspond to the same quantities at the begining of the white dwarf phase (after the H-reignition close to the knee in the HR-diagram). fThis 
sequence did not undergo a DIN, see text for details. 
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Table 2. Outburst properties of the DINs studied in this work. Unless stated otherwise, all sequences come from ZAMS progenitors with Z = 0.0001. Timescales are defined as follows; g 

cooling time at the moment of the CNO-flash (i~i), expansion time from the maximum energy release to the giant stage at logT c ff = 3.9 (tj), duration of the cool (logT fj < 4) giant ~ 

stage (ts) and contraction time needed to reach the pre-outburst luminosity (T4). fThese remnants were obtained by reducing the mass lost during the final AGB thermal pulse from 

that predicted by standard AGB wind prescriptions. fThis remnant was obtained by applying an artificially high wind during the fifth thermal pulse. i These DINs were obtained form g 

the O.53946M sequence by including overshooting (OV) in the convective zone generated during the CNO-flash. * In this case the intershell compositio n of the remnant was modified g. 

to resemble the C- and O- rich surface abundances of PG1159 stars. IM86 indicates the sequence computed by Iben & MacDonald l llben fc MacDonald 198^). 3 
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0.73 


6.6 x 10 


-3 


0.038 




1.1 x 10- 3 


O.85M 


O.53946M (w/OV, CO-rich) 


0.23 


0.62 


0.033 




0.081 




0.033 


IM86 (Z = 0.001) 


O.6M 


0.29 


0.71 


6.6 x 10 


-4 


7.5 x 10" 


-4 


8.1 x 10" 6 



Table 3. Outburst abundances of the DIN sequences presented in 



Tabled 



